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Abstract 

We discuss the necessary conditions for a two-level system in the presence 
of an isotropic band edge to be transparent to a probe laser held. The two- 
level atom is transparent whenever it is coupled to a reservoir constituted 
of two parts - a hat and a non-hat density of modes representing a PEG 
structure. A proposal on the reconstruction of the band edge prohle from 
the experimentally measured susceptibility is also presented. 
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1. Introduction 


The investigation of optical properties of atoms coupled to dissipative 
environments with a structured density of modes has been a topic of ac¬ 
tive research over the years [ill8|. Of particular interest is the discussion of 
atoms (impurities) embedded in two or three-dimensional periodic dielectric 
structures, known as photonic crystals j9-32|, since they allow control over 
the electromagnetic density of modes and the spatial modulation of narrow- 
linewidth (high-(5) modes, in both microwave and optical regimes j^. When 
these structures are used to create one or several forbidden frequency bands 
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they allow control or complete suppression of spontaneous emission, as well 
as absorption from those embedded impurities H, 1^, [l6l4l8l. 0 Ell. [23l-l32j| . 
It was particularly relevant the early observation that a two-level atom em¬ 
bedded in a PBG [Hi 12, 113,1 [l3 | could retain some population in the upper 
level, even when the transition frequency was in the transmitting band, being 
the final state a dressed state of the atom with a localized field mode, which 
lies in the forbidden band. More recently the attention has been shifted to 
quantum dots embedded in photonic crystals where each individual quan¬ 


tum dot can be seen as an “artihcial atom” (33l-35|. The important feature 
in any of those situations above is that the “atom” placed in such a structure 
interacts with the held modes in the propagating frequency band and in the 
forbidden photonic band gap (PBG) as well, giving rise to many interesting 
coherent phenomena such as the possibility of controlling non-markovian de¬ 
cay Bl E^ . localization of superradiance (l7l |. quantum interference effects in 


spontaneous emission |20, |28|, transpar ency to a probe held [21|, and squeez¬ 


ing in the in-phase quadrature spectra |25l |. 

The majority of contributions regarding radiative properties consider only 
spontaneous emission of two, three, four and hve level atoms embedded in 

SSlaJ, with only a few exceptions treat- 


a PBG structure 20, 22, 31 


ing absorptive and dispersive properties [21|, [3^. As an important example 


of this last case, the absorption and dispersion properties of a A-ty pe atom 


decaying spontaneously near the edge of a PBG was studied [2l|]. It was 


pointed out, within an isotropic PBG model, that the atom can become 
transparent to a probe laser held, even when other dissipative channels are 
present, suggesting that many surprising ehects in the absorption and dis¬ 
persion of atoms embedded in such structures can appear. Most of those 
ehects were considered inside model systems composed by three or more lev- 


36l. l37l |. while they were not proved to be 


els 0, 0 0 H 0 0, 0, [32 

strictly necessary. 

Pursuing this line we revisit the problem of transparency of an atom 
placed near an isotropic band edge |2l|, but consider the minimal situation 
of transitions between two-levels only. We show that for it to be transparent 
to a weak driven held, the two-level atom must be coupled to a reservoir 
constituted of two parts - a hat and a non-hat density of modes representing 
a PBG structure. Transparency is therefore an inner property of the reservoir 
engineering. As a side result of this approach we consider the related inverse 
problem considered in [13,00 on the possibility to obtain information 
about the band edge prohle from two-level temporal decay in such structure. 
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Here we show that is also possible to reconstruct the band edge characteristics 
directly from the experimentally measured susceptibility. 

This paper is organized as follows. In Section ([2]), we present the model 
considered and its stationary solution. In Section (|3]), the linear susceptibility 
is evaluated, and two models of isotropic band gap structures are analyzed. 
In Section dl]), is discussed how to reconstruct the band edge characteristics 
from the experimentally measured susceptibility. Finally, in Section ([S]) we 
conclude the paper. 

2. Model 

The system considered here is a two-level atom with excited and ground 
state |1) and |0), respectively and with transition frequency cjq. The atom 
is probed by a weak electric held with frequency cj detuned from wq by 5 = 
cu— cuq. The decay of the excited state is due to a coupling with vacuum modes 
described by a collection of harmonic oscillators with frequencies Um- In the 
rotating wave approximation and in the interaction picture the Hamiltonian 
of the system is given by 

H = |0) (1| + H.c.) + |0) (1| + H.c) , (1) 

m 

where = —pioF'o is the Rabi frequency, pio is the atomic electric dipole 
moment, and Eo is electric held amplitude. The Qm represents the coupling 
between the atom and the vacuum modes and and bm are the creation 
and annihilation operators for excitations in the reservoir, with m = A, k 
indicating a photon state with polarization A and momentum k. For sake 
of simplicity we assume and Qm as real. In the period of time t the state 
of total system, atom + reservoir modes, can be written as a superposition 
given by 

iV'(i)) = ao{t) |0, {0}) -h |1, {0}) + ^ |0, {m}). (2) 

m 

The coefficients ao(t) and ai(t) are the probability amplitudes to hnd the 
atom in the ground and excited state and the photon reservoir in the vacuum 
state, respectively, while the coefficient a^it) gives the probability amplitude 
to hnd the atom in the ground state and a single photon in the state m of 
the vacuum modes. Substituting Eq. ([2]) into the Schrodinger equation 
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containing the hamiltonian ([T]) and projecting into each state at the right- 
hand-side of Eq. ([2]) gives the following equations of motion for the time 
dependent coefficients ai and am, 


iao{t) = Qai{t), (3) 

iai{t) = Qao{t) - 6ai{t) + '^gme~"^‘^'"~‘^°~^^*am{t), (4) 

m 

iam{t) = ai{t). (5) 

Integrating Eq. ([S]) and eliminating the vacuum amplitude in the equations 
for ao(t) and ai(f) it follows that 


iao(t) = flaiit), 


idi{t) = ^lao{t) — 6ai{t) — i / K{t — t')ai{t')dt\ 


where the kernel, K{t — t'), is given by 


( 6 ) 

(7) 


K{t -1') 


E 2 

ym^ 

m 


( 8 ) 


All the information about the reservoir is contained in the kernel above, 
which will be dependent on the frequency distribution of the vacuum modes. 
As we are interested in effects strongly dependent on the reservoir modes 
distribution we keep the integro-differential equation ([7]) without any ap¬ 
proximation. The procedure to solve Eqs. ([6]) and ([7]) is straightforward by 
Laplace transform, which for the initial conditions ao(0) = 1 and ai(0) = 0 
results in 


ro(s) 


s — -I- G{s) 

s[s — i6 + G(s)] -|- ’ 


(9) 




—iVt 

s[s — i6 + G(s)] -|- ’ 


( 10 ) 


where Vo(s), Ei(s) and G(s) are the Laplace transforms of ao{t), ai{t) and 
— respectively. We must only assume two hypothesis to proceed with 
our calculations. The hrst one is about the reservoir - its memory function 
must allow the atomic system to reach a steady state, since we are interested 
in equilibrium properties at this regime. Secondly, we are interested in the 
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linear atomic response to the external field, so the coupling between field and 
atom is considered to be weak in such a way that ao{t) 1 for all times. 
This means that we are considering a perturbative solution for ai{t), which 
is linear in the driven field amplitude. This approximation does not affects 
considerably the results we are discussing, since the only effect of considering 
a first order term is to neglect power broadening and saturation effects due 


the laser field intensity In first order in Eq. flTU]) simplifies to 


Ei(s) = 


n 


s[is + (5 + fG(s)] ’ 


and the steady state solution is obtained by the limit procedure 
a.((^oo) = limsV,W = j-p^. 


( 11 ) 


( 12 ) 


where we assumed that the memory function is in such a way that the limit 
s —)■ 0 for G(s) shall exists. 


3. Susceptibility 

The induced polarization due the applied external field is given given by 

poo 

P{t)= (13) 

Jo 

where x{t) is the complex susceptibility, whose imaginary and real part are 
related to the atom absorption and dispersion of energy from the laser field, 
respectively j4l|. For a harmonic field E{t) = Eoe~'^^^ + h.c. the polarization 
becomes 

P{t) = 2Re {e-^‘^\{uj)Eo) , (14) 

where x{^) is the fourier transform of x{t). On the other hand, the atom 
polarization is obtained as an average of the atomic dipole moment, 

P{t) = (/i(t)) = 2Re (poi ao{t)al{t)) ^ 2Re (/ioiat(t)) • (15) 

From Eqs. (ITTll and (lT5l) . and considering N atoms by unit of volume, we 
obtain for the susceptibility at the stationary state, the following expression 

X(5) = (16) 
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which contains the information about the environment through the function 
G(0). Here, we assume that the atomic decay, i.e., the emission rate into 
the vacuum modes can be divided into two parts Q - the response 

due to the flat modes of vacuum plus a response due to structured modes. 
Following this assumption, we can write the kernel as 


K{t -1') = ^i(t -1') + k(t -1'). 


( 17 ) 


The first term in Eq. ffT7|) corresponds to a Markovian evolution, due to the 
coupling to a flat density of modes where 7 = 4a;g|/ioiP/3c^ is the Wigner- 
Weisskopf decay rate [i^ . The second term is the non-Markovian counterpart 
of the evolution, corresponding to the coupling with the structured density of 
modes imposed by the photonic crystal. The linear susceptibility now writes 


as 




6-q-iG*{oy 


(18) 


where G(0) is the limit t ^ 00 for the Laplace transform of the non-flat part 
of the kernel, K{t — t'). 

As a first example, we consider that the non-flat part is due to the coupling 
to an isotropic photonic band gap where the effective mass dispersion relation 
is = Ug + A(|k| — |ko|)^, with A a;g/|kop Uj, ll2|, |23|. In this case the 
kernel is given by 


d^/^P“d77/4+(<59-<5)h-t')] 

Kit - t') = ^ , -, t > t', 

TTtVti) 


(19) 


with = 2a;y^|/ioiP/3c^ and 6g = Ug — ujo- This model corresponds to a 
density of modes given by p{u)') = 0{u)' — ujg)f7rqu' — Ug, where 9{u' 
is the Heaviside function. The non-flat kernel then gives us 


OJn 


G{s) = 


^3/2g-m/4 

qs + i{ 6 g - 6 )' 


( 20 ) 


and the linear susceptibility reads 


x(^) ~ 


(s - ij/2) 


( 21 ) 
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Figure 1; Absorption and dispersion as a function of the detuning for the isotropic 
effective mass model. Solid line: Absorption; Dashed line: Dispersion. We set 6g = 2, 
/? = 1, and 7 = 1. The parameters are dimensionless. 

From Eq. (I2T]) it can be seen that the susceptibility is zero at 6 = 6g (or 
CO = Ug), and the atom is transparent to the probe laser held. To illustrate, 
in Fig. [T]we plot Re{x) and —Im{x), which corresponds to dispersion and 
absorption, respectively, as function of S, setting 6g = 2, (3 = 1, and 7 = 1. 
Besides transparency (zero absorption) at S = Sg = 2, it is also observed a 
strong deviation from the typical two-level absorption-dispersion curves in 
the Markovian approximation. Thus, the same transparency phenomenon of 
a three-level system in a A conhguration in the presence of a band edge 
is also possible for a two-level atom. 

In addition, is worthwhile to consider two extreme cases: {i) when the 
hat part of the density of modes vanish, and {ii) when the step-like function 
describing the band-gap changes in a smooth fashion, as in a real photonic 
crystal. In the hrst case, for S < Sg, except by a delta function absorption 
spike at the shifted atomic resonance frequency, the susceptibility given by 
Eq. fl^T]) has no imaginary part, and therefore, there is no absorption at 
these frequencies. For the second case we consider the one band model with 
a smooth density p{uj') = — Ug9{u' — Ug)/Tr{u' — Ug + e), where e is a 
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Figure 2: Absorption and dispersion as a function of the detuning around S = dg for the 
smoothed one-band isotropic model. Solid line: e = 1; Dashed line: e = 0.1; Dotted line: 
e = 0.01. We set = 2, /3 = 1, and 7 = 1 . The parameters are dimensionless. 


smooth parameter to avoid the singularity at u' = Ug B0- The kernel for 
this density is 

^ ( 22 ) 


and the susceptibility, 


i^/e + -f- i{6g - 5) ’ 


x(^) 




{S - i-f/2) {y/6g - S + ^/e) + /33/2 


(23) 


In this case the susceptibility is null at 6 = 6g in the strict case of e = 0. 
Otherwise, the absorption and the dispertion curves present a dip at S = 6g 
dependent on the value of e. In Fig. [2]we plot Re{x) and —Im{x) as function 
of 6 around ^ = Sg considering several values of e. The dip in the line shapes 
approach zero as the values of e goes to zero. 

As a second model for the reservoir we assume that the memory kernel 
has a non-markovian term due a two-band isotropic effective mass model for 









the photonic band gap structure. This model corresponds to a density of 
modes given by 


and resulting in 


_ 9{uja - u') 9{uj' - Ub) 

^ 27l\/Ua — Oj’ 2 tI^UJ' — CJfe’ 


(24) 


G(s) 


^3/2gm/4 

2 \/s + i{5a — 


where 6a = ooa — (^o, and 6b = oob — ^o- 
following expression 


+ 


^3/2g-i^/4 


i) 2^s + i(St-6)’ 


(26) 


Now the susceptibility is given by the 


^_ y( 6 „ - - 6 ) _^ 

(6 - n/2) V(6. - 5)(s, - i) + (-»y(4 - 6) + V(6„ - 6)) ' 

(26) 

In this case the susceptibility is null at 6 = 6a or 6 = 6b, and the atom 
becomes transparent at two different frequencies. The behavior of dispersion- 
absorption curves is illustrated in Fig. [3]where we plot Re{x) and —Im{x) as 
function of 6 and setting 5a = 1 and 6b = 2, (3 = 1, and 7 = 1 . Transparency 
at two different frequencies is observed. 


4. Band-edge profile reconstruction 

Since we have developed all the necessary ingredients for understanding 
the effects of the structured reservoir on the atomic transparency to the 
probe, we would like to discuss on a rather important related problem, which 
is the inverse problem of determining the characteristics of the band gap 
and its profile from experimental data. As pointed out by Nabiev 0 , the 
band gap prohle by can be determined from the experimental data on the 
temporal behavior of the atomic spontaneous decay. Indeed, this can also 
be done by a stationary measure through the susceptibility function. Under 
the conventional continuum limit for the reservoir mode distribution, and 
for long times, the Laplace transform for the non-flat contribution can be 
written as 

poo poo 

hmG'(s) = lim / duj'T{uj') / dr 

s-i-0 Jq Jq 

= 7iT{cj)-zV (27) 

Jo 
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Figure 3: Absorption and dispersion as a function of the detuning for the two-band 
isotropic effective mass model. Solid line: Absorption; Dashed line: Dispersion. We set 
(5o = 1 and = 2, /3 = 1, and 7 = 1 . The parameters are dimensionless. 


where r(a;') = g'^{u')p{u') represents the product of the coupling with the 
density of states for the non-flat sector, and we have used the identity ji^ 



7i6(u' — u — is) — iV -, 

UJ' — UJ — IS 


(28) 


and V means principal value. 

For the band edge prohle reconstruction it must be noticed that the hrst 
term in Eq. fl27)) is exactly taken in the external held frequency. Thus, 
when the external held frequency is varied, the reservoir frequency is probed, 
he., with the variation of the probe held frequency, in fact the band gap 
frequency distribution is scanned. If one assumes the general expression (ITHll 
and inverts G(0) as a function of the measured susceptibility, namely 


G(0) = 


7 

2 



A*(w) / ’ 


(29) 


and the band gap prohle can be reconstructed. 
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Remark that the present approach can be applied as well to a broad range 
of distinct situations, such as in the recent hndings on the coupling of atoms 
trapped in the near held of nanoscale photonic crystal cavities In this 
situation the present approach would be useful for probing the cavity density 
of modes through a susceptibility measurement. 

5. Conclusion 

In conclusion, we saw that it is possible to obtain transparency to a laser 
probe held on a two-level atom if it is considered the atomic coupling with 
a reservoir constituted by hat and non-hat densities of modes existing in a 
PEG. We have considered two isotropic band-gap models and analized the 
linear response to a weak optical held through the susceptibility function. We 
have also discussed the possibility of band edge prohle reconstruction via the 
susceptibility function knowledge. This can be an alternative to the method 
of reconstruction of the band edge prohle involving the measurement of tem¬ 
poral quantities as suggested by |l^ and implemented with single quantum 
dots embedded in a photonic crystal (See also Ref. ji^ for another 
method of reconstruction). In contrast the band edge prohle reconstruction 
here described can be realized in a steady state situation. Those hndings 
are particularly relevant for the emerging held of quantum nanophotonics 
[isj l. as well as in the investigation of nonlinear features with actual atoms 
trapped in the near held of nanoscale photonic crystal cavities [i^. The 
measured susceptibility can therefore be a valuable tool for band edge prohle 
reconstruction in those systems. 

This work was partially supported by CNPq and FAPESP through the 
Institute Nacional de Ciencia e Tecnologia em Informagao Quantica (INCT- 
IQ) and through the Research Center in Optics and Photonics (CePOF). 
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